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Stannary 

> 

Efficient  reanalysis  models,  which  provide  high 
quality  explicit  approximations  for  the  structural 
behavior,  are  introduced.  The  presented  algorithms 
are  based  on  a series  expansion  which  is  shown  to  be 
equivalent  to  a simple  iteration  procedure^  To  pre- 
serve efficiency,  only  methods  which  do *6t  involve 
matrix  inversion  have  been  considered.  yCmly  the  decom- 
posed stiffness  matrix,  known  from  exact  analysis  of 
the  initial  design,  is  required  to  obtain  the  approxi- 
mate expressions.  Two  approaches  of  accelerated  con- 
vergence are  proposed  to  improve  the  quality  of  the 
approximations : 

a)  An  approach  where  a- scalar  multiplier,  used  for 
scaling  of  the  initial  design,  is  chosen  prior  to 
the  solution  as  the  a elerating  parameter. 

b)  An  approach  where  information  gathered  during 
calculations  of  the  series  coefficients  is  used 
to  improve  the  convergence  rate. 

Numerical  examples  illustrate  the  efficiency  and 
the  quality  of  the  proposed  approximations.  A special 
attention  is  focused  on  reanalyses  along  a line,  a 
problem  typical  to  many  optimal  design  procedures.  The 
computational  effort  in  this  case  is  considerably 
reduced,  since  only  a single  independent  variable  is 
irvnlvpH  . 


1.  Introduction 

In  most  optimal  design  procedures  the  behavior  of 
the  structure  must  be  evaluated  many  times  for  succes- 
sive modifications  in  the  design  variables.  This 
operation,  which  involves  much  computational  effort, 
is  one  of  the  main  obstacles  i applying  optimization 
methods  to  large  structural  .terns.  Reanalysis 
methods,  intended  to  analyze  efficiently  new  designs 
using  information  obtained  from  previous  ones,  can 
broadly  be  classified  as  [1]: 

(a)  Direct  methods,  giving  exact  solutions  and  appli- 
cable to  situations  where  a relatively  small  pro- 
portion of  the  structure  is  modified  (for  example, 
only  a small  number  of  elements  are  changed). 

(b)  Iterative  methods  [2,31,  suitable  for  cases  of 
relatively  small  cham  i in  the  structure.  The 
known  solution  of  a gi.eti  design  is  usually  used 
as  an  initial  value  for  the  iterative  process. 
Problems  of  slow  convergence  rate  or  even  diver- 
gence may  arise  for  Urge  changes  in  the  design. 

(c)  Approximate  methods  [4-8] , usually  based  on  series 
expansion  and  require  less  computational  effort. 
One  problem  often  encountered  is  that  the  accuracy 
of  the  solution  may  not  be  sufficient.  Under 
certain  assumptions,  some  approximate  methods 

are  shown  to  be  equivalent  to  iterative  procedures. 

In  this  study  reanalysis  methods  for  optimum 
structural  design,  based  on  explicit  approximations  of 
the  structural  behavior  in  terms  of  the  independent 
design  variables,  are  presented.  Once  the  explicit 
model  has  beer,  introduced,  it  can  be  used  for  multiple 
reanalyses  of  designs  obtained  by  successive  changes 
in  the  variables.  The  presented  algorithms  are  based 
on  a series  of  expansion  which  is  shown  to  be  equiva- 


lent to  a simple  iteration  procedure. 

Two  conflicting  factors  should  be  considered  in 
choosing  an  approximate  behavior  model  for  a specific 
optimal  design  problem: 

(a)  The  computational  effort  involved,  or  the  effic- 
iency of  the  method. 

(h)  The  accuracy  of  the  calculations,  or  the  quality 
of  the  approximation. 

To  preserve  the  efficiency,  the  presentation  is  limited 
to  methods  which  do  not  involve  matrix  inversion.  Only 
the  available  decomposed  stiffness  matrix,  known  from 
exact  analysis  of  the  initial  design, is  required  to 
obtain  the  explicit  expressions.  However,  since  the 
proposed  models  are  based  on  a single  exact  analysis, 
the  accuracy  of  the  approximations  might  be  sufficient 
only  for  a limited  region. 

Two  approaches  of  accelerated  convergence  are 
proposed,  to  improve  the  quality  of  the  approximations: 

(a)  An  approach  where  a scalar  multiplier,  used  for 
scaling  of  the  initial  design,  is  chosen  prior 
to  the  solution  as  the  accelerating  parameter. 
Several  algorithms  for  selecting  the  value  of 
this  multiplier  are  proposed  and  their  merit  is 
demonstrated. 

fc)  An  approach  where  the  accelerated  parameters  are 
calculated  from  results  obtained  during  the 
solution  process.  Information  gathered  during 
calculation  of  the  series  coefficients  is  used 
to  improve  the  convergence  rate. 

Some  numerical  examples  illustrate  applications 
of  the  proposed  procedures.  A special  attention  is 
focused  on  reanalyses  along  a given  line  in  the  design 
space,  a problem  common  to  many  optimal  design  pro- 
cedui  s.  The  efficiency  and  the  quality  of  the  pro- 
posed approximations  are  demonstrated. 

2 . Problem  Statement 

The  displacement  analysis  equations  for  a given  design 
variables  vector  {X}  are 

[K]  {rf  = [R>  U) 

where  [K]  = stiffness  matrix  cori esponding  to  the 
design  {?J;  { R}  = load  vector  whose  elements  are 

assumed  to  be  independent  of  the  design  variables; 
and  tr}  = nodal  displacements  computed  at  { X}  . The 
elements  of  the  stiffness  matrix  [ K ) are  some 
functions  of  the  design  variables  {X}.  Assuming  a 
change  LA 5t>  in  the  design  variables  so  that  the 
modified  design  is 

{>5  = {*}  ♦ <AXl  (2) 

the  corresponding  stiffness  matrix  is  given  by 

[V]  = [K]  ♦ [AK]  (3) 

where  [A K j = the  matrix  of  changes  in  the  stiffness 
matrix  due  to  the  cl  uge  {Ax}  . 

The  object  in  this  study  is  to  present  explicit  models 
for  efficient  calculation  of  the  displacements  (r  } 
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corresponding  to  designs  { X>  , obtained  by  changing 
the  value  of  the  design  variables.  It  is  assumed  that 
the  displacements  i r } are  known  from  analysis  of  the 
initial  design.  Also,  [fc]  is  given  in  the  decomposed 
form 

Kl  = [Ul'lUl  (4) 

* 

where  [U]  is  an  upper  triangular  matrix. 

Approximations  along  a given  line  in  the  design  space 
are  often  required  in  optimal  design  procedures.  This 
problem  is  common  to  many  mathematical  programming 
methods  such  as  flasible  directions  or  penalty  function. 
A set  of  lines  (.or  direction  vectors  ) in  the  design 
space  are  determined  successively  by  the  optimization 
method  used.  In  each  of  the  given  directions  it  is 
usually  necessary  to  evaluate  the  constraint  functions, 
or  to  repeat  the  analysis,  many  times.  A line  in  the 
design  space  can  be  defined  in  terms  of  a single 
independent  variable  Y by 

(X}  = (X>  ♦ YlAXi  (5) 

* * 
where  { X}  is  the  given  initial  design,  {AX}  is  a given 
direction  vector  in  the  design  space,  and  the  variable 
Y determines  the  step  size.  Approximations  along  a 
line  require  much  less  computations  since  only  a single 
variable  is  involved. 

In  general,  the  elements  of  the  stiffness  matrix  are 
some  functions  of  Y.  dme  common  case  is  that  the 
modified  s;  ffness  matrix  can  be  expressed  as 

[X]  --  (K1  - f(Y)  [AK]  (6} 

In  truss  structures  where  (X)  are  the  cross-sectional 
areas  or  in  beam  elements  where  the  moments  of  inertia 
a'e  chosen  us  design  variables,  the  elements  of  the 
stiffness  matrix  are  linear  functions  of  V anu  bq.  (o) 
becomes 

IK]  * [KJ  > Y(AK]  (71 

If  the  elements  of  [K]  arc  functions  of  aZ^  ( Z i being 

the  naturally  chosen  design  variables  and  a.b  are  given 
cci stants)  we  may  use  the  transformation 

XA  = aZ*jJ  (8) 

and  obtain  the  linear  relationship  (7). The  expiession  of 
Hq_(g)is  suitably  for  example,  for  standard  joists  [9J  . 
Tn  cases  where  such  transformations  are  not  possible 
i lor  example,  in  frame  elements  where  the  stiffness 
matrix  is  a function  of  both  moments  of  inertia  and 
cross-sectional  areas),  still  linear  approximations 
may  be  used  for  the  nonlinear  terms  of  the  stiffness 
matrix  [10] . 


3 ■ Explicit  Behavior  Modi  s 
** 

The  analysis  equations  at  (x)  are 


m 

{r}  - [Rj 

(S) 

Based  on  Eq.  (3) , 

([K] 

+ 

t* 

7** 

J: 

ii 

70 

(1«) 

Premultiplying  by 

* - 1 

IK]"1  and  substituting 

* 

{rj 

1 

11 

(ID 

* 

* - 1 * 

IB] 

= [K]  1 [AK] 

(12) 

yields 

an 

* 1 8 j ) (V)  . {£} 

(13) 

Premu.lt ip'  •'ing  by 

([!)  + IB])  * and  expanding 

UI]  * [»] )' 1 - m - 18}  ♦ |B|  [B]1  + 

. (14) 

tq.  (13)  becomes 

<r}  = (flj  - [»]  + |B]2  - pV  a ...)(i*}  (16) 

Ihe  coefficients  of  this  series  can  readily  be  cal- 
culated. Defining 

* * k 

(rj)  - {I  r (16) 

f r2 } = - [ B] { r2 } (171 

etc.,  the  series  of  Eq.  (i5)  becomes 

* * * + 

{•■}  = lr>  <■  { r j > * {r^l  ♦ ...  (igj 

lor  the  given  triangularizat ion  of  bq.  (4),  the  cal- 
culation of  the  coefficient  vectors  { r } , (r^} , . . , 

requires  only  forward  and  back  substitutions.  The  cal- 
culation of  {?,},  for  example,  is  carried  out  as 
follows.  Substituting  tq,  (17)  into  bq.  (lo)  and  re- 
arranging gives 

IKJlrj)  * - [AX] t r } = (Rj}  (19) 

* 

he  first  solve  for  (P)  by  a forward  substitution 

[U]r{P}  = {R^  (70) 

* 

(r  j } is  then  calculated  by  the*  backward  substitution 
(UJiTj}  = {P}  I21j 

The  coefficient  vectors  {rj,{r  } etc.,  can  be  cal- 
culated in  a similar  manner. 

It  is  instructive  to  note  that  the  series  of  bq.  (15) 
is  equivalent  to  the  simple  iteration  procedure  [7,8) 

. (Ij  - [Bjtf1-”}  (22) 

where  k denotes  the  iteration  cycle  and 

(V(o)  } = (r  } (25) 

In  the  case  of  approximations  along  the  line  defined 
by  bq.  (5),  the  expression  of  Eq . (IS)  will  become 
explicit  function  of  {rl  in  terms  of  Y.  Assuming  the 
relationship  (6) , we  obtain 

{r}  « Ui]-f(Y)[l)  * f2(Y)[B|2  - ...){?}  (24) 

If  the  linear  dependence  of  Eq.  (7)  holds,  this  expli- 
cit expression  beconu  ; 

(r)  = ([I]-[B]Y  ♦ [Bj2Y2-[BJ  Vw  ...){*}  (75) 

or  (see  Eq.  (18)) 

{r}  = (r)+  { r x > Y * (r2}Y2  ....  (2b) 

This  equation  can  readily  be  used  for  mul  iple  re- 
analyses along  a line.  Also,  it  can  be  shown  that 
bq.  (26)  and  lay  lor  scries  expansion  of  the  displace- 
ments are  equivalent  [7,8].  Other  approximate  methods 
can  be  used  [7,8],  however,  these  usually  involve 
matrix  inversion. 

Klii le  the  methods  discussed  so  far  are  based  on  a 
single  exact  analysis  at  {X},  it  should  be  recognized 
that  better  approx -mat ions  could  be  obtained  if 
results  of  two  exact  analyses  (at  and  (10)  were 
considered.  Assuming  for  example,  quadratic  and 
cubic  interpolations,  respectively,  we  find 
* * 

{rl  = { rl  + (~}Y  - ({*r}  - {r}  - (||l)Y2  (27) 
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tr)  ~ [r) 


(Iy)Y 


1. 3 { V}  -3{r)  -2{|l}  - 

* * * 

* [2{rf  - 2{V>  ^ ) * C|y>)y3  (28) 

I'he  derivatives  can  readily  be  computed  by  several 

methods  [1],  One  possibility  is  to  differentiate  Eq.fl) 
with  respect  to  Y,  The  result  is 


f.K]  {|fj  - - [ff]<ri 


(29) 


in  which  both  {r},  and  [K]  in  the  decomposed  form  of 
Hq.  (4),  are  known  from  the  analysis.  Thus,  solution 

for  {—■}  involves  only  calculation  of  the  right  hand 

side  vector  of  Eq.  (29)  ami  forward  and  backward  sub- 
stitutions . 


t displacements  at  IX  1 is 

r<tV)  y 
r((X  })  V 


(39) 


This  result  indicates  that  the  error  in  the  approxima- 

tions  at  iXci  depends  only  on  the  ratio  r/r  (and 

not  on  & , the  distance  between  <x*l  and  (X  }).  The 

c 

combination  of  approximate  behavior  models  and  scaling 
can  be  used  to  introduce  efficient  optimal  design  pro- 
cedures 1 10 1 . It  will  be  shown  in  the  next  section 
how  selection  of  the  scaling  multiplier  a may  improve 
the  approximate  behavior  models. 

5 . Improved  Approximations  by  Scaling 


4 . Behavior  of  scaled  Designs 


Scaling  of  the  initial  design  (X)  to  obtain  a modified 

■30) 


design  {X  } is  giver,  by 
<x 


{X  } - a{X} 
a 

where  a is  a positive  scalar  multiplier,  if  the 
elements  of  the  stiffness  matrix  arc  assumed  to  be 
linear  functions  of  the  design  variables  then 

[Ka]  = a[KJ  (31) 

and  the  displacements  of  the  modified  design  are 
itq.  i 1 ) j 

(r  I - - (r)  (32) 

a u 

The  significance  of  this  relation  is  that  a given 
design  {£}  can  easily  be  scaled  by  modifying  & so 
that  any  desired  displacement  be  equal  to  a predeter- 
mined value, (an  operation  called  staling  of  the  design. 
(The  line  aix)  is  called  a design  line). 

In  optimal  design  problems  it  is  often  necessary^o 
find  the  design  (Xc)  along  the  design  line  ot{X}  , 
with  a displacement  r equal  to  its  limiting  value 


That  is  (see  Fig.  1) 
1 


(33) 

(34) 


In  cases  where  r is  calculated  by  an  approximate 
behavior  model  (such  as  Eq.  (15j),  we  may  evaluate  the 
accuracy  of  the  displacements  of  the  approximated 
design  (Xet  as  follows.  {X  ] is  givtn  by 

(X  } = aix] 


where  3 is  determined  from 


1 '* 
wr  r 
a 


(35) 

(36) 

3 = (37) 

ru 

tew  1 ** 

in  which  r is  the  approximated  value  of  r at  the 
point  (v ; . The  approximated  displacement  of  {)[  ) is 
r(!H.})  - ru  and  the  ex  >1  displacement  at  this 
point  is  (l;q.  (37)) 

rC(  JCJ  ) - r*  » rU  (38) 


The  modified  design  (X)  can  be  expressed  as  (see 
FiR.  2) 

{*X>  = (X)  * (AX)  = {X  } + (AX  ) (40) 

u a 

and  the  corresponding  stiffness  matrices  are  (Eq.(31)) 

[*K  1 = [ K J *[AK]  - [KJtlAKJ  = a i K ] + ) (41) 

Substituting  Eq.  (41)  into  Eq.  (9),  premultiplying  by 
[K]~I  and  rearranging  yields 


([I]  *x  [K)_1[AK  1)o(V}  i;> 


(42) 


Substituting  I A Kq 1 from  Eq.  (41)  into  Eq.  (42)  gives 


(II) 

Defining 


~ U)  [K]_1lAk])  afrl  - ir}  (43) 


. 1-q 
"o'  ~ a 


,-l. 


l-q 

a 


substituting  into  Eq.  (43)  and  expanding  ([!]+[»  ])  , 
w_  obtain  the  following  series  for  (V)  in  tcrmsJof  Qt 
(see  Eq.  (IS)) 


(r)  - i (!l]-(Ba]w[Baj2-|Ba]3v  ...){r; 


(45) 


For  a = l we  find  [Bu,  = [F 1 and  the  series  of  Eqs. 
(4S)  and  (15)  become  equivale.-t.  Different  dnection 
vectors  (AXa)  may  be  selected  in  the  plane  of  (K)  and 
{^0  for  various  a value.,.  While  it  is  usually 
difficult  to  predict  which  o wil'  irovide  improved 
convergence,  some  possibilities  e summarized  in 
Tabic  1.  In  cases  a,b,c,  the  value  of  a is  chosen  so 

that  the  resulting  direction  {AX  } is  perpendicular 
* **  a 
to  (X),  (X)  , and  to  the  bisector  of  angle  0 , 
respectively.  The  criterion  in  cases  d,e  is  chosen 
such  that  the  elements  on  the  principal  diagonal  of 
(AK  j or  the  second  term  in  the  series  of  Eq.  (45)  be 
equal  zero.  In  both  cases  the  multipliers  cu  aTe 
chosen  separately  for  each  displacement.  Results 
obtained  for  different  a v lues  will  be  compared  in 
the  numerical  examples  of  .ect.ion  7. 

The  scaling  multiplier  a effects  both  the  direction  vec- 
toT  {AX  } and  the  step  size  | AX  j,  where  (see  Eqs. 

(30)  and0 (40) ) 


(AX, ,1  = <V}  - a{X.) 


(4b) 


The  smallest  step  size  is  determined  by  (cose 
Table  1) 


That  is,  '.he  ratio  between  ihe  approximated  and  the 


a(x)  {AXul  = 0 


(47) 


7-19 


or,  after  rearranging 

» ■ 4#  <«) 

(X]T(x] 

and  the  corresponding  step  size  is 

|AXj  = l|xV  - a2|X|i)‘4  (49) 

Another  parameter  which  represents  the  step  size  for  a 
given  direction  is  the  angle  6,  where 

r-n*  ft  * VAT  rrn. 


large  k values. 

D^mamic  acceleration 

In  cases  of  slow  convergence  rale  dynamic  acceleration 
methods,  which  make  use  of  previous  terns  in  the  series, 
car.  be  employed.  If  the  iterative  process  has  a slow 
convergence  rate,  successive  errors  will  generally 
exhibit  an  exponential  decay  in  the  later  stages  of  the 
iteration  (see  Fig.  3).  Aitken's  S'  process  [13] 
is  one  approach  to  predict  the  asymptotic  limit  to 
which  the  predictions  for  each  displacement  r,  is 
tending.  Assume  the  extrapolation  expression  J 


evidently,  for  any  given  direction  (AX  } a better  con- 
vergence will  be  obtained  for  smaller  a 8 values. 

It  is  instructive  to  note  that  a scalar  multiplier  8 
can  be  chosen  instead  of  a such  that  (Fig. 2) 

(Xg>  = B{X]  (511 

Defining 

IBg]  2 (8-1) [I]  r BLBl  (52) 

we  nay  obtain  the  series 

(V)  * 8(H)  - F Bg J * [Bg]2  - (Bg)3  <-...)[£}  (S3) 

Comparing  Eqs.  (S2) , (53)  with  Eqs.  (44) , (45) , it  can 
be  observed  chat  identical  results  would  be  obtained 
by  bath  series  if  a = 1/8.  In  this  case  the  directions 
(AX(,}  and  (AXg}  are  parallel,  and  the  convergence 
rate  will  be  identical  for  any  given  6. 

6 • Convergence  Considerations 

Problems  of  slow  convergence  or  divergence  may  be  en- 
countered in  applying  the  series  of  Eq.  (15).  The 
series  converges  if,  and  only  if  [11], 

lim  [Bk]  = [0]  (54) 

k-»<*> 

A sufficient  criterion  for  the  convergence  of  the 
series  is  that 

IIBII  <_  1 (55) 

* * 

where  IIBII  is  the  norm  of  (B] . It  can  be  shown  that 

P(IBJ)  < IIBII  (56) 

* * 
i which  p( [B] ) is  the  spectral  radius  of  matrix  [B], 
uefined  as  the  lurgest  eigenvalue  |A,| 

D([B])  « (A  | = max  |A  | (57) 

i 

From  Eqs.  (55)  and  (56)  we  have 

P([B])<  1 (58) 

It  should  be  noted  that  the  existence  of  a norm  such 
that  IIBil  > 1 doe*,  not  preclude  the  convergence  of  the 


Some  procedures  have  been  proposed  to  predict  the 
eigenvalue  .\ j , One  possibility,  based  on  the  use  of 
Rayleigh  quotient,  is  [12] 


* T * . * * T » 

t'V  trk}  [rkrl} 

A.  “V  j J “ * „ v 

(rk)T  (rR>  Uy  (rk] 

* 

where  (r.  } are  th«  vectors  of  the  series  (soe  hq. 
(18)).  A*better  estimation  would  be  obtained  for 


where  a,  b,  c are  constants  and  k is  the  iteration 
number  (or  number  of  terms  in  the  series) . The  final 
solution  (k-«»)  is  determined  from  the  three  successive 
estimates 


(51) 

r(k> 

J 

= a + 

be'kc 

j 

* ft  ♦ 

(52) 

(k+2) 
r . 

- a ♦ 

be-^25c 

The  result  if 

(S3) 

#* 

r^ 

r(k*2)  . 

r . - a 
J 


or,  alternatively 


V.  - a . r!W)(S  (r<k+^  -r<k+1>) 

J J J 1 2 

r(X+l)  _ (k+2J 


If  Aitken's  acceleration  is  applied  at  the  wrong  time 
the  denominator  of  Eq,  (62)  could  be  zero  or  very 
small.  In  such  circumstances  the  method  will  either 
fail  to  yield  a prediction  or  else  give  a predicted 
value  which  is  grossly  in  error.  The  "WTong  time"  may 
be  interpreted  a*,  either  too  soon,  before  an  exponen- 
tial decay  is  established,  or  too  late  .(hen  rounding 
errors  affect  the  predictions. 

Jenr.ings[12]  proposed  a modified  version  cf  Aitken's 
acceleration,  used  by  several  authors  [14,15],  While 
the  prediction  of  Eq.  (62)  is  calculated  for  each  dis- 
placement r-  separately,  a common  acceleration  para- 
meter is  introduced  for  all  variables  in  the  modified 
method.  The  result  is 

•<*»•*>)  (65) 
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The  convergence  rate  is  govern  3d  by  the  magnitude  of 
Aj.  This  method  is  particularly  effective  when  onl.. 
one  eigenvalue  of  matrix  [S]  has  modulus  close  to 
unity.  It  is  po:  ible  to  apply  the  acceleration  after 
two  or  more  iterations  and  to  repeat  the  procedure 
frequently. 


7.  Numerical  Examples 

T»n  bar  truss.  The  truss  sho-n  in  Fig.  4 is  subjected 
to  a single  loading  condition  (all  dimensions  are  in 
kips  and  inches)  and  the  initial  design  is  (X)  = (6.0). 
The  following  three  cases  of  changes  in  the  design  were 
solved: 

Case  1:  {AX}T  - {4. 0,4. 0,4. 0,4. 0,6. 0,6.0, 

1.2, 1.2,1. 2, 1.2} 

Case  2:  {AX}T  = {-3.0 , - J.O ,-3.0,- 3.0 ,9.0 .9. 0,9.0 , 

9.0. 9.0.9.01 

Case  3:  {AX)T  * {18.0,18.0,18.0,18.0,18.0,18.0, 

36.0. 36.0.36.0.36.0} 

The  angles  0 for  the  three  cases  (Eq.  (50)) are  11.3°, 
30°,  and  11.5°,  respectively. 

Results  obtained  for  cases  1,2,  by  Aitken's  5*  method 
(Eq.  (62))  and  the  modified  acceleration  method  (Eq. 
(65)),  assuming  k = 2,3,4,  are  given  in  Table  2. 

While  the  approximations  for  case  1 are  excellent,  some 
errors  can  be  observed  in  case  2.  The  iteration  his- 
tory for  the  latter  case,  with  Aitken's  6*  method 
applied  after  iterations  4 and  6,  is  shown  in  Figi.5 
and  6.  It  can  be  seen  that  no  convergence  of  the 
vertical  displacements  could  be  achieved  without 
Aitken's  method.  Applying  scaling  by  the  five  methods 
of  Table  1 combined  with  Aitken's  process  (k  = 2,3,4) 
for  case  2 may  improve  the  convergence,  as  shown  in 
Table  3.  The  best  resu  s have  been  obtained  by 
methods  b and  d.  The  ei  ect  of  a on  the  spectral 
radius  (Eq.  (59))  is  illustrated  in  Fig.  7.  The  diver- 
ence  for  a=l  is  explained  by  the  relatively  large 
A j | value  ( | A ^ j =1 . S) . Assuming  0=1.96,  the  value  of 

| X ( is  reduced  to  0.75, 

The  effect  of  scaling  on  the  convergence  is  demonstra- 
ted in  case  3 (Table  4).  Applying  the  modified 
acceleration  method  (k=2,3,4),  no  convergence  could 
be  obtained  for  a=l.  Assuming  a=4.94  (method  c. 

Table  1),  the  convergence  is  fast;  a solution  very 
close  to  the  exact  one  is  obtained  for  three  terms  in 
the  series.  Similar  results  could  be  reached  with  the 
criterions  of  methods  a,b  in  Table  1. 

Forty-seven-bar  truss  ■ The  truss  shown  in  F:.g.  8 was 
solved  for  the  {following  data  (all  dimensions  are  in 
kips  and  inches):* 

Initial  design  {X}  = {0.5}. 


AX.  = 

0.5 

(i=l, ... ,8) 

AXi  = 

0 

(1=9,10,27-30,37-40,45-47; 

A3L  = 

0.2 

(i=ll-20, 41-44) 

AX.  = 

0,4 

(i*21-26) 

AX,  - 

-0.15 

(i=31-3b) 

Notation  of  the  displacements  is  as  follows:  rj,r,».- • 

are  the  horizontal  displacements  and  ^.r^,.-.  are  the 

vertical  ones.  Results  obtained  by  Aitken's  method 
and  the  modified  acceleration  method  (in  both  cases 
k=2,3,4)  are  given  in  Table  S.  It  can  be  noted  that 
despite  the  different  order  of  magnitude  of  the 
various  displacements,  relatively  small  errors  have 
been  obtained.  To  illustrate  the  effect  of  the  step 
size  Y,  the  displacements  rjg  r^ , r4J  have  been 

calculated  alcrg  the  line  {X}  * {X}>Y{Ai}  for 
Y = 0,2S,  0.50,  0.75,  1.0,  Results  obtained  after  our 
iterations  (k»4)  are  shown  in  Fig.  9.  The  effect  of 
Aitken's  methed  on  the  iteration  history  for  Y=1.0 
is  illustrated  ir.  Fig.  10  and  tne  evaluation  of  jn,j 
is  demonstrated  in  Fig.  11.  The  final  vn'uo  ]Aj|=  099 
explains  the  slow  convergence  rate  of  the  series. 


8.  Concluding  Remarks 

Approximate  behavior  models  for  efficient  reanalysis 
have  been  presented.  The  algorithms  are  based  on  a 
series  expansion  which  is  equivalent  to  a simple  itera- 
tion procedure.  A single  exact  a-.-lysis  is  sufficient 
to  introduce  the  series  coefficients  and  matrix  inver- 
sion is  not  required  throughout.  The  proposed  approxi- 
mations might  be  sufficient  only  for  a limited  region, 
in  the  neighborhood  of  the  initial  design. 

Two  approaches  have  been  proposed  to  improve  the  quality 
of  the  approximations : 

a)  An  approach  wh  re  a scalar  multiplier,  used  for 
scaling  of  the  initial  design,  is  chosen  prior  to 
the  solution  as  the  accelerating  parameter. 

Several  algorithms  for  selecting  the  value  of 
this  multiplier  are  proposed  and  their  potential 
for  improving  the  series  convergence  is  demonst- 
rated. It  is  shown  how  the  scaling  multiplier 
affects  both  the  direction  vector  in  the  design 
space  and  the  step  size. 

b)  An  approach  where  the  accelerated  parameters  are 
determined  from  results  obtained  during  the  solu- 
tion process.  Information  gathered  during  cal- 
culations of  the  series  coefficients  is  used  to 
introduce  extrapolation  expressions.  The  two 
methods  of  Aitken's  6J  process  and  a modified 
method  of  Aitken's  acceleration  are  presented. 

It  is  shown  how  these  methods  provide  high  quality 
results  in  cases  of  poor  convergence  rate  or  diver- 
gence cf  the  series. 

In  the  typical  problem  of  reanalyses  along  a given  line 
in  the  design  space,  the  methods  discussed  in  this 
study  involve  much  less  computational  effort.  Multiple 
reanalyses  along  a line  can  efficiently  be  introduced 
in  terms  of  a single  independent  variable. 
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Table  1:  Various  Possibilities  for  the  Selection  of  a 


Case  Criterion  for  Determining  a 


{AX^}  perpendicular  to  a{X} 


{AX^}  perpcndicuK, r to  {x} 


(AX  1 perpendicular  to  the 
bisector  of  0 


AK. . - 0 

ixa 


second  term  in  the  series 


Condition 

a 

a{X}l{AX  } = 0 

U)T{7} 

{x}r{x} 

{X}r{AX  } = 0 

m 

Ok 

txiYxi 

|ctX|  , |*X| 

/xV(x*^ 

v » r~V  ) 

** 

K 

AK.  . - V. . -a.K.  . = 0 

a.  = Jk 

rra  n i J ' 

1 K^i 

{b. } ‘ {r} 

(1-a.)  r.  + {B.}  {r}=  0 
11  1 

Uj  *■  1 + J 

r. 

i 

Table  2:  Results  (*100),  Ten-Bar  Truss. 


Rn 

53 

S3 

0. 

0. 

0.4 

57. 

1 0. 

0.813  0.310 
0.818  I 0.295 
0.876  0.310 


73  0.219 
73  0.220 
73  0 .219 


243  0.294 
231  0.287 
243  0.313 


k » 2,3,4. 


Table  3:  Results  (*100) , Ten-Bar  Truss  Case  2 Eq.  [62),  for  Various  a Values. 


3.041 

0.878 

0.915 

0.877 

0.756 

0.309 

0,309 

0.309 

0.310 

0.310 

0.876 

■ 

0.310 

0.194 

0.304 

0.335 

0.243 

0.303 

0.243 

0.222 

0.243 

0.303 

0.243 

k « 2,3,4 
See  Table  1, 
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